Meta-QTL analysis and identification of candidate genes for multiple-traits associated with spot blotch resistance in bread wheat

In bread wheat, a literature search gave 228 QTLs for six traits, including resistance against spot blotch and the following five other related traits: (i) stay green; (ii) flag leaf senescence; (iii) green leaf area duration; (iv) green leaf area of the main stem; and (v) black point resistance. These QTLs were used for metaQTL (MQTL) analysis. For this purpose, a consensus map with 72,788 markers was prepared; 69 of the above 228 QTLs, which were suitable for MQTL analysis, were projected on the consensus map. This exercise resulted in the identification of 16 meta-QTLs (MQTLs) located on 11 chromosomes, with the PVE ranging from 5.4% (MQTL7) to 21.8% (MQTL5), and the confidence intervals ranging from 1.5 to 20.7 cM (except five MQTLs with a range of 36.1–57.8 cM). The number of QTLs associated with individual MQTLs ranged from a maximum of 17 in MQTL3 to 8 each in MQTL5 and MQTL8 and 5 each in MQTL7 and MQTL14. The 16 MQTLs, included 12 multi-trait MQTLs; one of the MQTL also overlapped a genomic region carrying the major spot blotch resistance gene Sb1. Of the total 16 MQTLs, 12 MQTLs were also validated through marker-trait associations that were available from earlier genome-wide association studies. The genomic regions associated with MQTLs were also used for the identification of candidate genes (CGs) and led to the identification of 516 CGs encoding 508 proteins; 411 of these proteins are known to be associated with resistance against several biotic stresses. In silico expression analysis of CGs using transcriptome data allowed the identification of 71 differentially expressed CGs, which were examined for further possible studies. The findings of the present study should facilitate fine-mapping and cloning of genes, enabling Marker Assisted Selection.

in warm and humid areas of countries including Bangladesh, Nepal, Brazil, India, and Zambia 2 .The disease also affects the quality, texture, and colour of harvested wheat grains 3 .
Resistance against individual diseases in all major crops, including wheat, has largely been treated as a qualitative trait.Therefore, R genes for resistance have been identified and utilized for conventional methods of plant breeding.Spot blotch is no exception to this trend, so that four major genes, namely Sb1 to Sb4 following a gene-for-gene (GFG) relationship, are now known and have been largely used for resistance breeding; more recently, the Tsn1-ToxA system for sensitivity to this disease following an inverse gene-for-gene (IGFG) relationship has also been discovered [4][5][6][7][8] .During the last two decades, however, resistance against most diseases, including spot blotch, has been treated as a quantitative trait for genetic studies.These relatively recent studies included linkage-based interval mapping for the identification of Quantitative Trait Loci (QTLs) and LD-based GWA (Linkage Disequilibrium-based Genome Wide Association) studies for the identification of MTAs.As a result, a large number of QTLs/MTAs have been identified for each individual disease, including spot blotch.However, relatively few individual QTLs or MTAs have been utilized successfully for MAS, leading to the development of resistant cultivars.This is attributed to the non-availability of robust QTLs/MTAs for diverse genetic backgrounds and environments.MQTL analysis is a technique that provides a solution to this problem, so that a large number of MQTL studies for disease resistance in wheat have recently been conducted, either for resistance against individual diseases or for multiple disease resistance (MDR).Following are some examples of individual diseases: stem rust 9 , fusarium head blight [10][11][12][13] , tan spot 14 , stripe rust 15,16 , leaf rust 17,18 .Similarly, examples of MQTL analysis for MDR include two of our own studies for resistance against five diseases in one study 19 and resistance against all three rusts in another study 20 .
For spot blotch in wheat, more than a dozen studies involving interval mapping have so far been conducted 2 .However, for the reasons mentioned above, only sparingly QTLs have been recommended and used for the transfer of spot blotch resistance using MAS 21 .The use of MAS for breeding resistant cultivars will, however, be facilitated if robust MQTLs with high phenotypic variation explained (PVE) and narrow confidence interval (CI) are available.Keeping this in view, the present study was planned, where QTLs not only for SBR, but also those for some other associated traits were utilized for the development of MQTLs.These other traits included the following traits: (i) stay green (SG); (ii) flag leaf senescence (FLS); (iii) green leaf area duration (GLAD); (iv) green leaf area of the main stem (GLMS); and (v) black point resistance (BPR).These other traits were found to be associated with spot blotch in several QTL interval studies [22][23][24][25][26][27][28] .The present study is the first study involving MQTL analysis for spot blotch resistance and associated traits for the development of robust markers to be used for MAS.Candidate genes (CGs) associated with genomic regions occupied by MQTLs were also identified for a better understanding of the genetic basis of spot blotch resistance in wheat.

Bibliography search
A literature search involving QTL interval mapping for spot blotch and related traits was undertaken, which resulted in the identification of 228 QTLs from 28 studies involving 24 mapping populations (Supplementary Table S1).These QTLs were located on 20 different wheat chromosomes (except 6B; Fig. 1a).The maximum number of QTLs was available on chromosome 1B (41) followed by 2A (27) and 5A (23); the minimum number of only one QTL was available on 5D (Fig. 1a, b, c, e).

Consensus map and QTL projection
An integrated consensus map was prepared using markers that were utilized for interval mapping in all earlier studies.This map carried 72,788 marker loci (mainly SNPs, SSRs, and DArT markers; SNP = Single Nucleotide Polymorphism; SSR = Simple Sequence Repeat; DArT = Diversity Array Technology) (Fig. 1d, Supplementary Table S2).The number of markers on individual chromosomes ranged from a minimum of 591 (4D) to a maximum of 7,126 (1A).The genetic lengths of individual chromosomes ranged from 163.2 cM (4B) to 441.5 cM (5B) (Fig. 2) The minimum length between adjacent markers on a chromosome was 6.9 cM on chromosome 7D, and the maximum length (27.9 cM) on chromosome 1A (Supplementary Table S3).Only 84 of the above 228 QTLs could be used for projection; the remaining 144 QTLs did not have the complete information required for projection.Out of 84 projected QTLs, only 69 were useful (Table 1); the remaining 15 QTLs were singleton (each based on a single QTL) and therefore, were not utilised in downstream analysis.Among the 69 projected QTLs, the highest number of 23 QTLs were available for the trait GLAD (33.3%), followed by 20 (29.0%) for SBR and 19 (27.5%) for GLMS.

MQTL analysis
MQTL analysis gave 16 MQTLs distributed on 11 wheat chromosomes (Table 1).The number of QTL involved in individual MQTL ranged from a maximum of 17 QTLs in MQTL3, 8 each in MQTL5 and MQTL8 and 5 each in MQTL7 and MQTL14; the remaining MQTLs were each based on < 5 QTLs (Table 1).The number of traits controlled by individual MQTL ranged from 4 traits (SBR, GLAD, GLMS, and BPR) involved for MQTL5 to only single traits (SBR) in each of five MQTLs (MQTL2, MQTL4, MQTL11, MQTL15 and MQTL16) (Table 1).Out of the six traits, SBR was controlled by 12 (75%) out of 16 MQTLs, while GLMS was controlled by each of 9 MQTLs (56.2%);SG and BPR were each controlled by only 2 MQTLs (12.5%).The distribution of MQTLs on three subgenomes also differed with 8 MQTLs on A sub-genome, 6 on B sub-genome and 2 on D sub-genome (Fig. 2).Similarly, chromosome-wise distribution of MQTLs ranged from a maximum of 2 each on five chromosomes, namely 1B, 2A, 5A, 5B and 6A and a minimum of 1 each on several other chromosomes (Fig. 2).Mean PVE % for individual MQTLs ranged from 5.4 to 21.8% (mean 13.5), the logarithm of odd (LOD) ranged from 3.2 to 11.5 (mean 5.3) and the CI for individual MQTLs ranged from minimum 1.5 cM (MQTL2) to a maximum on 57.8 cM (MQTL1) with a mean of 16.5 cM (Table 1).This amounted to 4.6 times reduction in the length of CI in original QTLs.The average length of CI on chromosomes 2A was reduced by 27.7 fold and that on chromosome 5A was reduced by 26.4 fold, followed by 12.7 and 8.4 fold reduction in MQTLs on chromosomes 1B and 4B.The spot blotch resistance gene, Sb1 4 was associated with MQTL16 (Table 1, Fig. 2).

Selection of major MQTLs
The above 16 MQTLs were also subjected to selection of the most important MQTLs using the following four criteria: (i) high frequency of QTLs (ranging from 4 to 17) associated with the corresponding individual MQTL; (ii) high PVE% (ranging from 10.7 to 21.8), (iii) high LOD score (ranging from 4.7 to 6.9) and (iv) relatively narrow CI (ranging from 2.4 cM to 10 cM).This exercise gave the following four major MQTLs: MQTL3, MQTL4, MQTL5, and MQTL8.These four major MQTLs were later also found to be associated with CGs that were particularly relevant to resistance against biotic stresses (Table 2).

Candidate genes
The genomic regions associated with the above 16 MQTLs carried 516 unique CGs, which encoded ~ 508 proteins; 411 of these proteins are already known to be associated with tolerance to a variety of biotic stresses, including resistance/sensitivity to spot blotch.(Supplementary Table S5).www.nature.com/scientificreports/Gene Ontology (GO) of CGs.The gene ontology (GO) of CGs suggested occurrence of GO terms belonging to all three widely known classes of functions, including cellular, molecular, and biological processes.The details of these functions are presented in Fig. 3a.
The fold change of the DECGs ranged from 5.4 fold upregulation to -5.0 fold downregulation.The maximum difference in expression of DECGs was observed in the pair of susceptible and resistant genotypes at 24 hpi and 48 hpi (R24 vs. S24 and R48 vs. S48).In these comparisons, four DECGs were significantly downregulated in four MQTLs (MQTL1, MQTL2, MQTL9, and MQTL15).Similarly, nine DECGs were upregulated in six MQTLs (MQTL3, MQTL4, MQTL5, MQTL6, MQTL12, and MQTL13).In both comparisons, sample taken from a resistant genotype was used as the control, suggesting that more genes are upregulated in S genotype relative  to those, which are upregulated in R genotypes.It also suggested that resistance is achieved by silencing many genes that are responsible for susceptibility.
An the above analysis involving comparison of R and S genotypes, DECGs also led to the identification of 58 unique domains in which 10 protein domains were associated with biotic stress resistance.Among the downregulated DECGs, the encoded proteins belonged to following four main categories: (i) Nucleic acidbinding, (ii) RING-type zinc-finger, (iii) Palmitoyltransferase, and (iv) GDP-L-galactose/GDP-D-glucose phosphorylase.Similarly, in upregulated DECGs, the following six protein domain were identified: (i) PPMtype phosphatase domain superfamily, (ii) DUS-like, FMN-binding domain, (iii) AOX protein family, (iv) Thioredoxin-like superfamily, (v) NB-LRR domain superfamily, and (vi) VQ domain (Fig. 3b, Supplementary Table S6).
In comparison of expression within the same genotype (either compatible interaction of incompatible interaction) at different intervals, some DECGs showed distinct expression patterns associated with compatible interactions (S0 vs. S24, S0 vs. S48, and S24 vs. S48), while others were known to be associated with incompatible interactions (R0 vs. R24, R0 vs. R48 and R24 vs. R48) (Supplementary Table S6).The functional analysis revealed distinct trends in protein domain functions between the two reaction types.In the compatible response, genes with downregulated expression were enriched for the following domains: (i) NAD(P) binding, (ii) carbon-nitrogen metabolism, (iii) specific regions of PRORP proteins, (iv) protein degradation (ubiquitinlike), and (v) protein synthesis (peptidyl-tRNA hydrolysis).In contrast, in the upregulated response involved in incompatible interactions, following domains were involved: (i) phosphoglycolate phosphatase-like domain, potentially linked to phosphate metabolism; (ii) splicing factor 3B subunit 5 and (iii) CYSTM protein domains.These domains suggested potential role of RNA processing and transmembrane signaling in defense response (Fig. 3b, Supplementary Table S6).

Discussion
In recent years, molecular breeding involving the use of DNA markers for MAS has become routine in conventional plant breeding.This is particularly true, when improvement of quantitative traits is involved.However, for MAS, robust markers that are suitable to be used for MAS are not always available.In such studies, MQTLs derived from QTLs identified in earlier studies may sometimes prove useful.For spot blotch and related traits, the present meta-QTL study is the first of its kind.The study utilized all known QTLs not only for spot blotch resistance but also those associated with five other related traits, which often accompany spot blotch resistance (for identity of five traits, see above).
As reported in results, in the present study, using 228 QTLs that were available from published literature, only 84.0 (36.8%) could be projected onto the consensus map, because only these had the complete information needed for MQTL analysis.This gave 31 MQTLs, of which 15 MQTLs were based each on a single QTL and were therefore rejected, thus leaving only 16 MQTLs derived from 69 QTLs.The remaining 144 QTLs could not be projected due to one of the following reasons: (i) there were no flanking markers that the original and consensus www.nature.com/scientificreports/maps shared, and (ii) the CIs were relatively large 19 .Notably, SG trait was the only trait besides spot blotch resistance that was directly targeted in QTL mapping within the context of spot blotch 23 .In previous studies [24][25][26][27][28] , the remaining traits (except BPR), namely GLAD, GLNS, FLS, were generally also subjected to mapping of tolerance against abiotic stress factors in wheat.A 14.2 -fold reduction in the number of genomic areas or QTLs linked to spot blotch resistance in wheat was achieved through the identification of 16 MQTLs from 228 QTLs.This projection success percentage is less than that of previous studies on MQTL analysis conducted for different wheat disease resistances, where the following projection rates were available: 44.0, 60.62, 66.6, and 75.2% 15,16,19 .The present study also adds to the examples of MT-MQTLs involving disease resistance, because 12 of the 16 MQTLs were MT-MQTLs.Earlier MQTL studies in wheat have been carried out for a variety of traits, including tolerance against abiotic and biotic stresses.For instance, Liu et al. 14 identified 20 MQTLs for tan spot using 106 QTLs, Amo et al. 17 reported 35 MQTLs for leaf rust using 128 QTLs, and Jan et al. 15 discovered 61 MQTLs for stripe rust using 184 QTLs.MQTLs for multiple disease resistance (MDR) were also identified in two of our own studies (see later for some details).
Although our focus in the present study was primarily on SBR, a five other traits known to be associated with SBR were also included.Among these associated traits,, SG trait was of special interest, since it has the unique characteristic of delaying senescence in leaves and other plant parts, enabling them to maintain high chlorophyll levels and sustain photosynthetic activity for an extended period.The association of SG with spot blotch has also been reported in wheat, where a significant negative correlation (0.73) between SG and AUDPC (a measure of spot blotch severity) has been reported 23,37 .In a number of studies, SG has often been used as a selection criterion for tolerance against abiotic stresses also, including heat and drought 16,[24][25][26]28,38 . Thus,this adaptive trait (SG) is also particularly useful in adverse environments, since it facilitates an improved grain-filling process 27 .The genotypes that are characteristics of SG and SBR also exhibit tolerance/resistance against heat/drought, suggesting a shared physiological mechanism for response to abiotic and biotic stresses including SBR 39 .
Since six different traits were used in the present study, it was also possible to identify MT-MQTL, which can be used for improvement of more than one traits using the same MQTL (see later in Discussion).In the published literature, only two other studies from our own work were available, which involved identification of MT-MQTLs, both involving resistance against a number of diseases.In one study, Saini et al. 19 reported MQTLs for MDR involving five unrelated diseases (STB, SNB, FHB, KB, and LS), while Pal et al. 20 reported MDR MQTLs for only the three rusts, including leaf rust, stem rust, and stripe rust.In the present study, out of 16 MQTLs, 12 involved more than one trait and were therefore described as MT-MQTL (involving 2-4 traits); no MT-MQTL was available for 5 or 6 traits (Table 1).Among these 12 MT-MQTLs, on the basis of priority, three MT-MQTL are being recommended for MAS involving SBR with other associated traits (see later).
The merit of MQTLs relative to the corresponding QTLs used for MQTL analysis, depends on four key attributes, including PVE, LOD score.CI and the number of QTLs involved in each individual MQTL.Following are some details of comparison of these four attributes in MQTL and QTLs used for MQTL analysis, suggesting that MQTL showed an improvement over QTLs: (i) The PVE of the MQTLs ranged from 1.6 to 55.2% (with a mean of 13.6%), while that of the QTLs ranged from 5 to 10% in 72 QTLs, followed by 10-15% in 60 QTLs, and 15 to 20% in 39 QTLs.This suggested that PVE% improved in MQTLs.(ii) The LOD scores of MQTLs ranged from 3.0 to 11.0; in contrast, in QTLs, the LOD score ranged from 2 to 4 in 83 QTLs, from 4 to 6 in 75 QTLs, and 6-8 in 48 QTLs, suggesting that the LOD score also improved in MQTL.(iii) The CI in MQTLs largely ranged from 1.5 to 20.7 cM (except five MQTLs, where it ranged from 36.1 to 57.8 cM) as against 0.62-272 cM for the QTL used for MQTL analysis, suggesting a major reduction in the length of the CI, thus making these MQTL relatively more robust.(iv)Similarly, 8 MQTLs were each based on two QTLs, and MQTL3 was based on 17 QTLs, suggesting that the number of QTLs associated with individual MQTLs is large (at least in some MQTLs), thus improving the credibility of the utility of MQTLs (Table 1).
The results of the present MQTL study on SB in wheat can also be compared with those available for several other wheat diseases, including three rusts 9,15,16,18 , fusarium head blight (FHB) 10 and tan spot 14 .Following are the some details about the MQTLs for disease resistance identified in other studies: (i) The PVE ranged from 4.7 to 51.0% for three rusts, 6.2-22.3%for FHB, and 6.3-27.0%for tan spot 14 .(ii) The LOD score for MQTL ranged from 3.0 to 28.0 for three rusts 20 , from 1.0 to 62.0 for FHB 12 , and from 3.9 to 17.0 for tan spot 14 .(iii) The CI in these earlier studies ranged from 0.04 to 83.5 cM for three rusts 20 , < 10.0-95.0cM for FHB 10 and 0.9-8.1 cM for tan spot 14 .For MT-MQTLs for MDR in earlier studies, the PVE ranged from 2.2% to 51.2%, the LOD score ranged from 2.9 to 48.4 and the CI ranged from 0.04 to 15.2 cM 19 .
For validation, a study of colocalization of MQTLs identified in the present study with MTAs identified in published literature was also undertaken.For this purpose, MTAs were available for only SBR, SG, and BPR; for other traits, no GWAS were available.For these three traits, 13 (> 81.0) of the 16 MQTs were colocalized with GWAS-MTAs (Supplementary Table S4.These results can be compared with similar validation studies undertaken earlier.In these earlier studies, co-localization ranged from 38.7 to 90.5% 17,20,[40][41][42][43][44] .The high rate of colocalization suggests the presence of causal polymorphisms in MQTLs identified in the present study. Another component of the present study was the identification of CGs associated with MQTLs.The distribution of CGs ranged from none in MQTL7 to a maximum of 157 in MQTL2, with an average of ~ 33 CGs per MQTL (Supplementary Table S5).This data suggests a potential role for these CGs in spot blotch resistance, warranting further investigation.The R domains in proteins encoded by CGs are known to be involved in plant defense mechanisms, strengthening the hypothesis of their association with disease resistance.In wheat, similar studies on CGs associated with MQTLs were earlier conducted for traits like fusarium head blight 12 , stripe rust tolerance 15,16 , leaf rust 17 , tan spot resistance 14 , MDR 19,20 .
Among a large number of CGs identified to be associated with MQTLs in the present study, 71 DECGs were available, which encode proteins carrying domains with known roles in plant defense.Some of these genes were downregulated (FC < − 2.0), while others were upregulated, when a resistant genotype was compared with www.nature.com/scientificreports/ a susceptible genotype at two different durations after inoculation (R24 vs. S24 and R48 vs. S48).The protein domains involved in downregulation included the following: nucleic acid-binding domain, RING-type zincfinger domain, and palmitoyltransferase domain.Similarly, the upregulated genes encoded proteins carrying the following domains: AOX protein family, AOX protein family, NBS-LRR and VQ protein.These results are also in agreement with the results of several earlier studies conducted either in Arabidopsis or in wheat [45][46][47][48][49][50][51][52] , suggesting that these DECGs deserve further detailed study.
In the present study, we also validated 13 (> 81.0%) of the 16 MQTLs using the previously conducted GWAS.This also allowed selection of the following for MQTLs based on high PVE, LOD value, and lower CIs: MQTL3, MQTL4, MQTL5, and MQTL8.Genomic selection models can also utilize these markers to enhance the accuracy of resistance prediction.On the basis of present study, we also recommend further research to clone and functionally characterize the identified candidate genes (CGs).These genes have the potential to be used to improve wheat resistance but require validation through techniques like gene cloning, reverse genetics, or omics approaches.
In the present study, among the four R genes (Sb1 to Sb4) known for spot blotch, Sb1 is the only Sb gene, which colocalize with one of the 16 MQTLs, namely MQTL16, suggesting that the list of MQTLs identified is certainly not exhaustive and exclusive and that there must be more SB MQTLs (associated with Sb2, Sb3 and Sb4) to be identified in future.
Among the 12 MT-MQTLs also, three MT-MQTLs (MQTL3, MQTL5 and MQTL8) are recommended for MAS on the basis of criteria outlined above for selection of MQTLs.MQTL16 carrying the R gene Sb1 is also recommended for MAS.Only future QTL and MQTL studies will allow identification of MQTLs associated with the remaining R genes, namely Sb1, Sb3 and Sb4, Some of the DECGs identified during the present study may also be utilized for further studies.

Summary
In the present meta-QTL study, 228 QTLs were available from earlier QTL interval mapping studies, involving spot blotch and five related traits.Only 69 QTL could be projected leading to the identification of 16 MQTLs; 12 of these MQTLs were MT-QTLs, each involving 2-4 traits, suggesting that the same meta-QTL may be used for improvement of more than one trait.Superiority of the MQTLs over corresponding QTLs was demonstrated on the basis of improved PVE%, LOD scores and reduced length of CI 11 .As many as 13 of the 16 MQTs were colocalized with GWAS-MTAs, thus placing s higher level of confidence in these MQTLs.A large number of CGs were also identified, one of them differentially expressed.The integration of desirable alleles from MQTL3, MQTL4, MQTL5, and MQTL8 can provide high and stable spot blotch resistance.MQTL3, with 17 initial QTLs for four traits including spot blotch disease resistance QTLs, may be particularly useful in developing spot blotchresistant cultivars through MAS.MQTL5 and MQTL8 each have eight initial QTLs.All five of these MQTLs meet the preferred criteria mentioned earlier and can be easily transferred into susceptible genotypes using MAS.Other MQTLs may also confer resistance against spot blotch disease after introgression into susceptible wheat backgrounds.

Literature survey for spot blotch related QTLs
For the identification of MQTLs, 228 QTLs were available from 21 different studies, with 24 mapping populations related to spot blotch resistance (SBR) and five other related traits.The related traits included SG, FLS, GLAD, GLMS, and BPR.GLAD and GLMS are both related traits, the former representing duration and the latter representing the extent of green leaf area; these have been treated as separate traits in published literature and therefore retained as two different traits in the present study also (Supplementary Table S1).

Development of consensus map
A consensus map was developed using 72,788 markers, including SNPs, DArT, SSRs, and EST-based markers, that were available in five distinct genetic linkage maps [53][54][55][56][57] .The software Lp merge was used for developing the consensus map 58 .Similar consensus maps were earlier developed and utilized to find the MQTLs for different traits in wheat 15,59 .The QTLs and the flanking markers for different traits were also used in the preparation of the consensus map.

QTL projection and MQTL analysis
Out of 228 QTLs, only 84 had all the information necessary for MQTL projection.These QTLs were used for projection, and meta-analysis was conducted using BioMercator V4.2 60 .When a marker's genetic position was uncertain, the QTL was projected on the consensus map, and the markers closest to the QTL were accepted as associated markers.Wherever the CI of an individual QTL was not available (there were 5 QTLs for which CIs were not available), the CI (95%) was estimated utilizing the following two equations, where N is the size of the population: (i) CI = 163 ÷ (N × R 2 ) for the QTLs identified using RIL populations 61 , (ii) CI = 530 ÷ (N × R 2 ) for F 2 and backcross (BC) populations 62 , and (iii) CI = 287 ÷ (N × R 2 ) for the QTLs detected using doubled haploid (DH) population 9 .
QTLs with inadequate information (with no PVE values, or no LOD scores or no information on genetic position, etc.) were eliminated.Following two different methods were used during the analysis, depending upon the available number of QTLs: (i) When the number of QTLs per linkage group was not more than 10, the method proposed by Goffinet and Gerber 63 was used; in this approach, we tested all possible combinations based on the Akaike Information Criterion (AIC) to determine the number of underlying meta-QTLs that best fit the available results, and (ii) when the number of projected QTLs per linkage group was > 10, the two-step analysis proposed by Veyrieras et al. 64 was employed; this is the best fit meta-QTL model based on the following parameters: (i) AIC, (ii) corrected AIC, (iii) corrected AIC with a penalty factor of 3, (iv) Bayesian Information Criterion (BIC), and (v) approximate weight of evidence criterion (AWE).The statistical algorithms and techniques furnished in the software are provided by Sosnowski et al. 65 .

Validation of MQTLs with GWAS-based MTAs
MQTLs identified in the present study were validated with MTAs reported from 21 independent GWA studies, including 11 GWA studies on SBR, 6 on SG and 4 on BPR.No GWAS studies were found for other traits used in the present study.These GWA studies utilized populations of hexaploid wheat only (e.g., spring and winter wheat), with population sizes ranging from 101 66 (BPR) to 1384 27 (SG), and phenotyped at one or more locations across five different countries (India, Mexico, Pakistan, USA for SBR, and China for BPR).The details of population size, associated diseases, genotyping platform, number of markers used, and MTAs detected in these studies are available in Table 3.The physical positions of each significant and stable marker associated with the trait were obtained from the respective studies or the JBrowse WHEAT URGI database (https:// urgi.versa illes.inra.fr/ jbrow seiwg sc/) and CerealsDB www.cerea lsdb.uk.net/ cerea lgeno mics/ Cerea lsDB/ index NEW.Table 3. Summary of genome-wide association studies (GWAS) in wheat on spot blotch (SBR), black point (BPR) and stay green (SG) used in the present study.

Figure 1 .
Figure 1.The data utilised in the MQTL analysis.(a) frequencies of QTLs used to identify MQTLs on each of the 20 wheat chromosomes (except 6B); (b) frequencies of QTLs with different LOD scores; (c) frequencies of QTLs for different traits; (d) density of molecular markers on consensus map involving 21 wheat chromosomes; (e) frequencies of QTLs with different PVE% values.

Figure 2 .
Figure 2. MQTLs are distributed on various wheat chromosomes, with MQTLs shown by blue and GWAS validated MQTLs indicated by red.

Figure 3 .
Figure 3. (a) Gene ontology (GO) for differentially expressed candidate genes (DECGs) identified in the present study; (b) Heat maps exhibiting 71 DECGs that encoded important proteins and are associated with the MQTLs.The scale shown on the top of left of the heat map shows the fold change value that varies in the heatmaps proteins encoded by DECGs. https://doi.org/10.1038/s41598-024-63924-w

Table 1 .
Summary of the results of MQTL analysis for resistance to spot blotch disease in wheat.*Sb1 gene is located within the interval of markers. MQTL-

Table 2 .
Candidate genes associated with each of four major MQTLs and their encoded proteins.